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Analysis of finite temperature phase transition using level spacing * 
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Let B be the largest spacing between adjacent eigenvalues of the Polyakov loop. We propose to employ the 
distribution of B as an order parameter for the finite temperature phase transition in SU (N) lattice gauge theories. 
Using smeared links to reduce ultraviolet fluctuations, we carry out a test for the gauge group SU{3). 



The lattice regulator is a powerful non- 
perturbative framework to study the ther- 
mal properties of gauge theories. Polyakov 
loops provide a useful order parameter for the 
study of the finite temperature confinement- 
deconfinement phase transition of non-Abelian 
SU{N) gauge theories. Polyakov loops are not 
properly renormalized quantities and in four di- 
mensions are strongly affected by ultra-violet fluc- 
tuations. It is well known that in the continuum 
limit at fixed temperature the expectation value 
of an unrenormalized Polyakov loop vanishes even 
in the deconfined phase, due to divergent self- 
energy contributions to the free energy. 

For iV < 3, a lattice theory with Wilson action 
undergoes a bulk cross-over, separating the lattice 
strong- and weak-coupling regimes [3] . The bulk 
cross-over is characterized by a steep increase in 
the average plaquette value and occurs at a fixed 
lattice coupling Pb- This steep increase occurs 
when the probability of 1 x 1 Wilson loops to 
have eigenvalues near —1 decreases dramatically 
(for = oo a true gap opens up) . For N > 5 the 
cross-over becomes a discontinuity and there is 
a phase transition, but no symmetry breaks and 
there is no continuum counterpart to this transi- 
tion. However, the bulk cross-over can affect the 
finite temperature transition [4] on coarse lattices. 

The Polyakov loop L is the parallel transporter 
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Figure 1. Eigenvalue distribution of the plaquette 
in the confined and deconfined phase. 

around a closed loop in the "temperature" direc- 
tion and gets multiplied by a Zj^ phase under cer- 
tain action preserving changes of link variables in 
the path integral. In the confined regime (the dis- 
ordered phase of a gauge theory) the eigenvalues 
of the Polyakov loops are randomly distributed 
over the unit circle. At some critical Pc the Zj^ 
symmetry gets broken and the eigenvalues of the 
Polyakov loop favor one of the N roots of unity 
on the unit circle. (3c scales with Nt, the length of 
the lattice in the temperature direction. In order 
to describe the physical finite temperature tran- 
sition it is necessary to have Pc > Pb- In our 
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SU{S) simulations the plaquette distribution al- 
ready had an apparent gap for Nf = A over the 
whole range of (3 that we studied as seen in Fig. 1 . 

Let the gauge invariant eigenvalues of a given 
Polyakov loop L be e**J [j = 1,2, ...,7V), 
-TT < 6*1 < 02 < ... < On < tt. 
The maximal level spacing B is defined as 
max{27r-6'Ar + 6ii, 9n - On-i, 6*2 -6*1}. 
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Figure 2. Effect of APE smearing the links. 

The distribution of B changes through the fi- 
nite temperature transition. It is neither affected 
by the particular location on the unit circle where 
the spectrum condenses, nor by tunneling be- 
tween the different Zn vacua. If the eigenval- 
ues 9j are randomly distributed on the unit circle 
B will be of order ^ 27r/iV; when the Polyakov 
loops get ordered, the concentration of angles at 
one location creates an effective gap, increasing 
B. Random matrix theory can be used to de- 
scribe the distribution of B in the confined phase 
where it has one peak value. This is no longer 
true in the deconfined phase where the distribu- 
tion tends to have two distinct peaks. The devel- 
opment of two peaks in the deconfined phase is 
signaled by the single peak distribution develop- 
ing a more populated tail as one goes through the 
phase transition. Deep in the confined phase one 
can compute the average value of B by averaging 
over a very large number of random matrices. In 
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Figure 3. {Re{L)) and (B) vs. (i 

the case of S'C/(3) we get {B.^^^^^^^) = 0.9564(1) 
and (Sj.j^q(Jqj-^) ^ for ^ oo in the case of 
SU{N). 

Following [1], we denote by Re(L) the pro- 
jection of the trace of L onto the nearest 
axis. The critical f3c is generally determined 
by locating the peak in xl = (Re(L)Re(L)) — 
(Re(L))(Re(I/)) . Projection onto the nearest Z^ 
works reasonably well for TV = 3 but can become a 
problem close to the transition for larger N . The 
observable B does not suffer from such a prob- 
lem since it is defined as the biggest level spacing 
without any reference to its location on the unit 
circle. 
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Susceptibility of the real part of the Polyakov loop 
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Figure 4. Susceptibilities of {Re{L)) and {B) . 

We performed various Monte Carlo simulations 
to study the finite temperature transition using 
the standard Wilson gauge action. A generic 
SU{N) gauge fields update was implemented, 
even though the present results are restricted to 
SU{'i). We used the heat-bath algorithm fol- 
lowing Creutz-Kennedy-Pendleton [5] , where the 
update of one link corresponds to updates by 
all the N{N - l)/2 SU{2) subgroups. In our 
Markov chain a step between two measurements 
consisted of one heat-bath sweep plus one over- 
relaxation sweep to avoid being trapped in a lo- 
cally metastable state. Each over-relaxation up- 
date was done on the full SU{N) group [6]. We 
considered lattices with Nt — 4: and L = 12, 16, 24 
and typically we collected at least 15,000 mea- 



surements. The errors were estimated with the 
jack-knife method. 

In order to suppress ultra-violet fluctuations 
we used APE smeared links following [7]. We 
smeared both the time-like and the space-like 
links on the lattice in parallel and repeated the 
process Nt times. Subsequently, we used the re- 
sulting links to compute the Polyakov loops. In 
Fig. (2) one can see that smeared links provide 
a largely enhanced signal in the Polyakov loop 
susceptibility. 

The biggest level spacing turned out to be a 
good order parameter showing transitions consis- 
tent with traces of Polyakov loops as one can 
see in Fig. (3) (the horizontal line represents 
^^random)' limiting value as /3 0). By 
studying xb in Fig. (4) we infer fic — 5.6925(25) 
and this is consistent with [1]. Both L = 16 and 
L = 24 show well defined peaks in xb and our 
estimate of Pc was done without resorting to re- 
weighting and by working at relatively small vol- 
umes. The lack of monotonicity seen in Fig. (4) 
in the deconfincd phase is due to the system be- 
ing trapped into one of the vacua for a long 
time during the measurement history. 
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